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Abstract 

The Gamow states describe the quasinormal modes of quantum systems. It is 
shown that the resonance amplitude associated with the Gamow states is given 
by the complex delta function. It is also shown that under the near-resonance 
approximation of neglecting the lower bound of the energy, such resonance am- 
plitude becomes the Breit-Wigner amplitude. This result establishes the precise 
connection between the Gamow states, Nakanishi's complex delta function and 
the Breit-Wigner amplitude. In addition, this result provides another theoretical 
basis for the phenomenological fact that the almost-Lorentzian peaks in cross 
sections are produced by intermediate, unstable particles. 
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1 Introduction 

Resonances appear in all areas of quantum physics, in both the relativistic and non- 
relativistic regimes. Resonances are intrinsic properties of a quantum system, and they 
describe the system's preferred ways of decaying. Experimentally, resonances appear as 
sharp peaks in the cross section that resemble the Breit-Wigner (Lorentzian) lineshape. 

The Gamow states are the natural wave functions of resonances, and they were 
introduced by Gamow in his paper on a-decay of radioactive nuclei [1]. Since then, 
they have been used by a number of authors, see e.g. [2-29]. Likewise the bound states, 
the Gamow states are properties of the Hamiltonian, and they are associated with the 
natural frequencies of the system. The usefulness of the Gamow states is attested 
by the remarkable success of the Gamow Shell Model [18, 22, 24-26, 28] and similar 
nuclear-structure formalisms [16,17]. 



Because resonances leave a quasi-Lorentzian fingerprint in the cross section, and 
because the Gamow states are the natural wave functions of resonances, the resonance 
amplitude associated with the Gamow states must be related to the Breit-Wigner am- 
plitude. The purpose of this paper is to show that the resonance amplitude associated 
with the Gamow states is proportional to the complex delta function, 5(E — zr), and 
that such amplitude can be approximated in the near-resonance region by the Breit- 
Wigner amplitude. More precisely, we will show that the transition amplitude from 
a resonance state of energy z R to a scattering state of energy E > 0, A(z& — > E), is 
given by 

A{z R —> E) = iV2^M K 5(E - z R ) ~ - ^L-= 1 , E > , (1.1) 

V 27T tj — ZR 



where 



A/^ = ires[S(£)W = ir R , (1.2) 



S denotes the S matrix, and r R denotes the residue of S at zr. In addition, we will 
see that the lower bound of the energy (threshold) is the reason why this amplitude is 
not exactly but only approximately given by the Breit-Wigner amplitude. 

Section [2] provides a quick summary of the most important properties of the Gamow 
states, along with some basic phenomenological properties of resonances. The proof 
of (11. ip is provided in Sec. [3j The conclusions are included in Sec. HJ 

For the sake of clarity, we shall prove Eq. (11.11) using the example of the spherical 
shell potential for zero angular momentum. However, as explained in Appendix [A], the 
result is valid for any partial wave and for spherically symmetric potentials that fall off 
faster than exponentials. Finally, in Appendix [Bj we provide a thorough characteriza- 
tion of the complex delta function and its associated functional, since they have rarely 
appeared in the literature. 



2 Basics of resonances and Gamow states 

Resonance peaks are characterized by the energy E R at which they occur and by their 
width Tr. The resonance peak is related to a pole of the S matrix at the complex 
number zr = En — iTR/2, because the theoretical expression of the cross section in 
terms of the S matrix fits the experimental cross section in the neighborhood of E R , 
see Eqs. O and ([£3]) below. 

When the peak is too narrow and its width cannot be measured, one measures the 
lifetime tr of the decaying particle. Decaying systems follow the exponential decay 
law, except for short- and long-term deviations. 

Although a decaying particle has a finite lifetime, it is otherwise assigned all the 
properties that are attributed to stable particles, like angular momentum, charge, spin 
and parity. For example, a radioactive nucleus has a finite lifetime, but otherwise it 
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possesses all the properties of stable nuclei; in fact, it is included in the periodic ta- 
ble of the elements along with the stable nuclei. Similarly, most elementary particles 
are unstable, and they are listed along with the stable ones in the Particle Data Ta- 
ble [30] and attributed values for the mass, spin and width (or lifetime). Thus, stable 
particles differ from unstable ones by the value of their width, which is zero in the 
case of stable particles and different from zero in the case of unstable ones. Hence, 
phenomenologically, unstable particles are not less fundamental than the stable ones. 

A priori, resonances and decaying particles are different entities. A resonance refers 
to the energy distribution of the outgoing particles in a scattering process, and it 
is characterized by its energy and width. A decaying state is described in a time- 
dependent setting by its energy and lifetime. Yet the difference is quantitative rather 
than qualitative, and both concepts are related by 

r R = -, (2.i) 

though in most systems one can measure either tr or Tr, but not both. 

Theoretically, however, the relation (12.11) is usually justified as an approximation, 
trTr ~ h, as a kind of time-energy uncertainty relation. For a long time, it was not 
possible to experimentally check whether the relation ( 12. II) is exact or approximate, 
since the lifetime and width could not be measured in the same system. This changed 
with the measurements of the width [31] and lifetime [32] of the 3p 2 P^/2 state of 
Na, which provide a firm experimental basis that Eq. (12.11) holds exactly, not just 
approximately. Thus, resonances and decaying systems are two sides of the same 
phenomenon. 

Although the resonance peaks in the cross section resemble the Lorentzian, the 
resonance lineshape does not coincide exactly with the Lorentzian. Two features of 
the cross section reveal so. First, the maximum of the resonance peak never occurs 
at E = En, whereas the maximum of the Lorentzian occurs exactly at E = E-&. And 
second, the Laurent expansion of the S matrix around the resonance pole, 

S(E) = -^- + B(E) 7 (2.2) 

& — Zr 

which produces the Lorentzian peak in the cross section [33], 

° ~ {e - E R y + (r R /2)2 ' (2 ' 3) 

is valid only in the vicinity of the resonance pole. Because ( 12.21) and ( 12. 3ft are valid only 
in the vicinity of the resonance energy, the Lorentzian lineshape is just a near- resonance 
approximation to the exact resonance lineshape. 

Because the Lorentzian does not coincide exactly with the resonance lineshape, the 
Breit-Wigner amplitude cannot coincide exactly with the resonance amplitude. One 
can reach the same conclusion by using the point of view of decaying states as follows. 
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The Breit-Wigner amplitude yields the exponential decay law only when it is defined 
over the whole of the energy real line (—00,00) rather than just over the scattering 
spectrum (see e.g. [34]). Because in quantum mechanics the scattering spectrum has a 
lower bound, the Breit-Wigner amplitude would yield the exponential decay law only if 
it was defined also at energies that do not belong to the scattering spectrum. Thus, the 
Breit-Wigner amplitude is incompatible with the exponential decay law, and therefore 
cannot coincide with the exact resonance/decay amplitude. 

Mathematically, the Gamow states are eigenvectors of the Hamiltonian with a com- 
plex eigenvalue z R = E R — iT R /2, 

H\z R ) = z R \z R ) , (2.4) 

and, in the radial position representation, they satisfy a "purely outgoing boundary 
condition" (POBC) at infinity: 



(r\z R ) ~ e V( 2ffl /Skr ) asr-^oo. (2.5) 

The time-independent Schrodinger equation ( 12 Ah subject to the POBC ( 12.51) is equiv- 
alent to the following integral equation of the Lippmann-Schwinger type: 

z R - H + 1O 

where Ho is the free Hamiltonian and V is the potential. Since Eq. (I2.6P also yields the 
bound states, the Gamow states are a natural generalization to resonances of the wave 
functions of bound states. The bound and resonance energies obtained by solving (I2.6P 
coincide with the poles of the S matrix. 

The time evolution of a Gamow state is given by 

e- iHt ' h \z R ) = e- iE ^ h e- r ^^\z R ) , (2.7) 

and therefore the Gamow states abide by the exponential decay law. Because the 
eigenvalue of Eq. (12 .4p is also a pole of the S matrix, Eq. (12. 7p implies that Eq. (12. ip 
holds. In this way, the Gamow states unify the concepts of resonance and decaying 
particle, and they provide a "particle status" for them. 

Furthermore, since one can obtain both the bound and the resonance energies from 
Eq. (12. 6p . or from the poles of the S matrix, resonances are qualitatively the same as 
bound states. The only difference is quantitative: The Gamow states have a non-zero 
width (i.e., finite lifetime), whereas the bound states have a zero width (i.e., infinite 
lifetime). 

An important feature of the Gamow states is that they form a basis that expands 
any wave packet p + , see e.g. review [35]. The basis formed by the Gamow states is not 
complete though, and one has to add an additional set of kets to complete the basis. 
In a system with several resonances, we have that 



= yy^kX*^) + / dEe- iEt / h \E+){ + E\cp + ) 



(2.8) 
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where z n = E n — iT n /2 denotes the nth resonance energy. In this equation, the sum 
contains the resonance contribution, whereas the integral contains the background. For 
simplicity, we have omitted the contribution from the bound states. The main virtue 
of resonance expansions is to isolate each resonance's contribution to the wave packet. 

Resonance expansions allow us to understand the deviations from exponential de- 
cay [36]. In the energy region where one resonance R is dominant, the expansion (12. 8p 
can be written as 

(p+(t) = e- iz * t/h \z R )(z R \<p + ) + background(R) , (2.9) 

where the term "background(R)" contains all contributions not associated with the res- 
onance R, including those from other resonances. Because "background(R)" will always 
be nonzero, there will always be deviations from exponential decay. The magnitude 
of these deviations depends on how well we tune the system around the resonance 
energy: The better we tune the system around the Gamow state \z R ), the smaller 
"background(R)" will be. Note that "background(R)" is the analog to the background 
B(E) of the expansion fl2T2l . 



3 Proof 

3.1 Preliminaries 

The proof of (II. ip presented below is a straightforward application of the theory of 
distributions. Rather than working in a general setting, we will use the example of the 
spherical shell potential, 

< r < a 

7(x) = V{r) = { V a<r <b (3.1) 
b < r < oo , 

and restrict ourselves to the s partial wave. 

In order to prove (11.11) . we need to recall that in quantum mechanics, the transition 
amplitude from one state to another is given by the scalar product of those states: 

A(z R ~^ E) = (- E\z R ) , E>0, (3.2) 

where (~E\ is the "out" bra solution of the Lippmann-Schwinger equation. 

For the potential (13. ip . the "out" Lippmann-Schwinger eigenfunction reads, in the 
radial position representation, as 

(r\E-) = x"(r; E) = N(E) ^| , E e [0, oo) , (3.3) 
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where N(E) is a delta-normalization factor 

N(E) = 

x(r; E) is the regular solution, 



1 2m/h 2 



vr ^/2m/h 2 E 



(3.4) 



sin(Ax) < r < a 

X (r; E) = { j!(E)e iQr + J 2 (E) e - iQr a<r <b 
J 3 (E)e ikr + J 4 {E)e~ ikr b < r < oo 



(3.5) 



the wave numbers k and Q are given by 

k = 



2m 



E , Q 



2m 



(E - V ) 



h 2 ' ^ V fr 2 

and J±(E) are the Jost functions, 

J+(E) = -2iJt(E) , J-(E) = 2i Js(E) . 



(3.6) 



(3.7) 



The resonance energies zr produced by the potential (13. ip coincide with the zeros of 
With each resonance energy zr, we associate a Gamow eigenfunction «(r; zr): 



(r\z R ) = u(r;z K ) 



where 



h) sin(fc R r) 



m 



< r < a 



Ml X ^ ^gi e iQnr + J(g e-iQar a < r < 6 ( 3 .8) 



^3(fcR) 



6 < r < oo 



R 



2m 2m 

zr , Qr = \/ — (zr - Vo) , 



and where A/r is given by Eq. (jl.2p . From Eqs. (13.31) and (13. 8p we obtain 

1 



X (r;z R 



iV^rA/R 



w(r; z R j 



(3.9) 



(3.10) 



As shown in [37], the analytic continuations of the Lippmann-Schwinger kets -and 
therefore also the Gamow kets- are well defined as antilinear functionals over the space 
of test functions ip~ for which the following quantities are finite: 



dr 



nr 



e ^ 2 /2(i + H) n 'i[>-(r) 



n, n' = 0,1,2, 



1 + nr 

The action of the Gamow ket |zr) on the test functions ip~ is explicitly given by 



(3.11) 



{fp \ z r) 



dr (ip \r) (r\z-&) 
6 



drip (r)*u(r; zr) . 



(3.12) 



For the sake of clarity, we need to introduce a special notation that will specify when 
we are working in the energy representation: Whenever we work in such representation, 
we will add a hat to the corresponding quantity. For example, the energy representation 
of a wave function ip~(r) will be denoted by ip~(E). 

By using the operator {/_ of [38], one can obtain the energy representation of ip~(r), 

POO 

i~(E) = (U_ip-)(E) = drif)-(r)x~(r;E)*. (3.13) 

Jo 

As shown in [37], when ip~{r) satisfies ( 13. lip , the analytic continuation of [ip~ (E)]* 
exists. We shall denote such analytic continuation by ip~(z*)*: 

POO 

$-( z *)*= dri)-(r)*x~{r]z) . (3.14) 
Jo 

Not onlv ip~{z*)* exist, all the test functions (z*)]* are analytic at the resonance 
energieso As explained in Appendix [Bj this means that the antilinear complex delta 
functional at the resonance energies z R can be defined as 

S ■ $ i > C 

r — » (rK)=^ R (r):=?-(4r, 

where $_ CX p is the space of test functions That is, 5 2R associates a test function 
with the value that the analytic continuation of takes at z R . 

3.2 The Gamow state and the complex delta function 

In order to obtain the equality of Eq. (11.11) . we are going first to denote the energy 
representation of the Gamow ket as 



%> = C/-kR>. (3.16) 



Then, it follows that 



$ 1%) = $ \ u -\ z r) 
= (U j J~\z R ) 



\zr) 

POO 

/ dr[il)-(r))*u(r;z R ) bv (I3T2D 

Jo 

POO 

i^M R dri/j-(r)* X "(r;z R ) by <^M> 

Jo 

iV2^ Af R $~ (z R )* by (J3HD 

iV2^ AT R (^-|5 2R ) by dSZSJ) . (3.17) 



lr The poles of [V> (z*)]* are located on the first sheet of the Riemann surface, that is, on the upper 
half of the wave-number plane. 
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This equation proves that the energy representation of the Gamow functional, \z R ), is 
proportional to the antilinear complex delta functional, \S ZR ). 

In the energy representation, the identity I can be written as 

dE\E-)(-E\=I, (3.18) 

where \E~) denotes the energy representation of \E~). By inserting (13.181) into the 
first and the last terms of (I3.17p . we obtain 

/»oo /*oo 

/ dE$-\E-)(-E\zj L )= / dEiV2^Af R $-\E-)(-E\5 ZK ). (3.19) 
Jo Jo 

Because Eq. (I3.19P holds for any ip~ , it follows that 

(-E\z~) = i^N R (-E\5 ZR ) = \V2^N K 5{E - z R ) , (3.20) 

which, after dropping the hat notation and using Eq. (13.21) . becomes the equality in 
Eq. (HID. 



3.3 The Gamow state and the Breit-Wigner amplitude 

In order to obtain the approximation of Eq. (II. ip . we are going to obtain first the 
transition amplitude A(z R — ► E) from a resonance state of energy z R to a scattering 
state of energy — oo < E < oo, 

A(z R -*E) = (-E\z R ), Ee (-00,00). (3.21) 

Even though a quantum system can only decay to a scattering state of energy E > 0, 
we are going to ask the system to pretend that it can also decay to negative energies. 
Mathematically, this is equivalent to ask the system to pretend that its scattering 
spectrum runs from —00 to +00. Physically, it is equivalent to ignore the effect of the 
lower bound of energy E = 0. Calculating (13.211) . i.e., forcing the system to decay also 
to negative energies, needs a regulator. The regulator we will use is e~ iaE , where a > 
and E has zero or negative imaginary part. The reason why we use this regulator is 
that for complex z, the wave functions ip~(z*)* grow slower than e' Im< -v^W^)l 2 j n ^he 
lower half plane of the second sheet [37]. More precisely, Proposition 3 in [37] shows 
that for each n = 1,2,... and for each (3 > 0, there is a C > such that in the lower 
half plane of the second sheet, 2p~(z*)* is bounded by 

^ I Im( v / 2?Ti/ft 2 z )| 2 

\r(zT\ < C N 1/4| 1 + Z |n e " • ( 3 - 22 ) 

This estimate implies that e~ iaz ip~(z*)* tends to zero in the infinite arc of the lower 
half of the second sheet, 

lim e- iza ^-(z*)* = , a > . (3.23) 
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In its turn, the limit (I3.23P enables us to apply Cauchy's theorem to obtain 



^ 1 r°° ^ i 

^-(4r = lim-— / dEe- iaE ^(ET——, (3.24) 

a^O 27T1 > /_ 00 E - Z R 

where the integral is performed infinitesimally below the real axis of the second sheet. 
By multiplying both sides of (I3.24p by i^/2nAf R , and by recalling (I3.17p . we obtain 

dEe- iaE ^(Ey(-l)^=- • (3.25) 

.oo V27T E - z n 

In the bra-ket notation, Eq. (I3.25P reads as 

/oo 
dEe- iaE $-\E-)(-E\z R ). (3.26) 
-oo 

Comparison of (13.251) with (I3.26P yields the following expression for the amplitude f)3.2ip : 
A(z R —> E) = - ^-= 1 , E e (-oo, oo) . (3.27) 

V 27T £j — Z R 



Thus, if the scattering spectrum was the whole real line, the resonance amplitude would 
be exactly the Breit-Wigner amplitude. However, because the scattering spectrum has 
a lower bound, the resonance amplitude is not exactly the Breit-Wigner amplitude. 
Only when we can neglect the effect of the threshold, the resonance amplitude coincides 
with the Breit-Wigner amplitude: 

A(z R - E) ~ A(z R -> E) , (3.28) 

which is the approximation on the right-hand side of (11.11) . In particular, when the 
threshold can be ignored, the complex delta function becomes for all intends and pur- 
poses the Breit-Wigner amplitude. 

It should be stressed that the amplitude (I3.2ip is not physical, because in (I3.2ip 
the energy E runs over the whole real line rather than over the scattering spectrum. 
However, such unphysical amplitude helps us understand what the physical amplitude 
-the complex delta function- is, by allowing us to see how the resonance would decay 
if the scattering spectrum was the whole real line. 



3.4 Further remarks 

Aside from phase space factors, cross sections are determined by the transition ampli- 
tude from an "in" to an "out" state, A(Ei — > Ef). If |E ± ) denote the "in" and "out" 
solutions of the Lippmann-Schwinger equation, then 

A(E t - E f ) = (-E f \E+) = S(E t )5(E f - E t ) . (3.29) 



9 



\A(, R ^E,)\^L^-— — — — ; (3.30) 



If we imagine now that instead of an initial state \Ef) we had an unstable particle 
|zr), the transition (decay) amplitude A(z R — > Ef) would be given by (11. ip . Using the 
approximate decay amplitude of (11.11) . one obtains the following approximate decay 
probability: 

1 -A/rI 2 i_ 

that is, the decay probability of a resonance is given by the Lorentzian when the 
effect of the threshold can be ignoredJl Because the almost-Lorentzian decay probabil- 
ity (I3.30P coincides with the almost-Lorentzian peaks in cross sections, resonances can 
be interpreted as intermediate, unstable particles. 

Finally, it is worthwhile to compare the Gamow states with the states introduced 
by Kapur and Peierls [39]. As mentioned above, the Gamow states are eigenfunctions 
of the Hamiltonian that satisfy the POBC (12. 5ft at infinity; the wave numbers involved 
in the POBC (12.51) are complex and proportional to the square root of the complex 
eigenenergies of the Gamow states; such complex eigenenergies are the same as the poles 
of the S matrix, and they do not depend on any external parameter or energy. By 
contrast, the Kapur- Peierls states are eigenfunctions of the Hamiltonian that satisfy a 
POBC at a finite radial distance r , where r is such that the potential vanishes for r > 
r ; the wave numbers involved in the POBC satisfied by the Kapur-Peierls states are 
real and proportional to the square root of the real energy of the incoming particle; the 
POBC satisfied by the Kapur-Peierls states makes them and their associated complex 
eigenenergies depend on tq and on the real energy of the incoming particle; also, the 
complex eigenenergies of the Kapur-Peierls states are not the same as the poles of the 
S matrix. Thus, the Kapur-Peierls states do not seem to be related to the standard 
Breit-Wigner amplitude, because such amplitude does not depend on r and its complex 
energy does not depend on the energy of the incoming particle. 



4 Conclusions 

Since resonances leave an almost-Lorentzian fingerprint in the cross section, and since 
the Gamow states are the wave functions of resonances, the decay amplitude provided 
by a Gamow state should be linked to the Breit-Wigner amplitude. In this paper, 
we have found that the precise link is given by Eq. (jl.ip . and we have interpreted 
this result by saying that the resonance amplitude associated with a Gamow state is 
exactly given by the complex delta function, and that the Breit-Wigner amplitude is 
an approximation to such resonance amplitude, which approximation is valid when we 
can neglect the effect of the threshold. Thus, Eq. (II. ip establishes the precise relation 
between the Gamow state, Nakanishi's complex delta function and the Breit-Wigner 
amplitude. In addition, Eq. (11.11) affords another theoretical argument in favor of 

2 A much more detailed study of the dependence of the cross section (and expectation values of 
observables) on the Breit-Wigner amplitude can be found in e.g. [2,8,14]. 
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interpreting the almost-Lorentzian peaks in cross sections as intermediate, unstable 
particles — resonances are real (as opposed to virtual) particles, in accordance with 
resonance phenomenology. 

As is well known, the actual resonance lineshape of cross sections can be very 
different from a quasi-Lorentzian one, due to the effect of thresholds, other resonances, 
or extra channels. The usefulness of (11. ip does not lie in predicting the exact shape 
of the cross section, but rather in identifying what contribution to the cross section 
comes from each pole of the S matrix. In particular, although the equality in Eq. ( 11. ip 
is always exact, for practical purposes the approximation in Eq. ( 11. ip is useful only for 
narrow resonances. 

When we add Eq. (II. ip to the other known properties of the Gamow states, we 
see that such states have all the necessary properties to describe resonance /unstable 
particles: 

• They are associated with poles of the S matrix. 

• They exhibit the correct phenomenological signatures of both resonances (almost- 
Lorentzian lineshape) and unstable particles (exponential decay), and they pro- 
vide a firm theoretical basis for ( 12.11) . 

• They are basis vectors that isolate each resonance's contribution to a wave packet. 

Acknowledgments 

The author thanks Casey Koeninger for encouragement, and Alfonso Mondragon for 
enlightening criticisms. The author thanks Alfonso Mondragon also for his precise 
explanation of the relation between the Gamow and the Kapur-Peierls states. 

A Generalizations 

Equation (ll.ip is not valid only for the spherical shell potential (13.11) but actually holds 
for a quite large class of potentials. The reason can be found in well-known results of 
scattering theory [40,41]. As explained in [41], page 191, partial wave analysis is valid 
whenever the spherically symmetric potential satisfies the following requirements: 

I. V(r) = 0(r~ 3 ~ e ) as r — > oo. 
II. V(r) = 0(r- 3 / 2+e ) as r -> 0. 

III. V(r) is continuous for < r < oo, except perhaps at a finite number of finite 
discontinuities. 
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These conditions are, however, not sufficient to guarantee that the S matrix S(E), 
the Jost functions J±(E) and the Lippmann-Schwinger eigenfunction x~( r ; E) can be 
analytically continued into the whole complex plane. Such analytic continuation is 
guaranteed when we replace condition I by the more stringent 

I. V(r) falls off faster than exponentials as r — > oo, 

as stated throughout Chapters 11 and 12 of [41], and in Chapter 5 of [40], especially 
in Theorem 5.3.2. Thus, when V(r) satisfies I-III, even though we may not know their 
exact analytic expressions, we know that S(E), J~±(E) and x~( r i E) can be analytically 
continued into the whole complex plane and that the Gamow eigenfunction u(r; %) 
is well defined. Moreover, since in the asymptotic region r — * oo the expressions of 
u(r; Zji) and x~( r ; E) for any potential satisfying I-III are the same as the expressions 
of u(r; zr) and x~{ r '->E) for the spherical shell potential in the region r > b (with 
different expressions for the Jost functions) , the general proof goes through exactly the 
same lines as the proof for the spherical shell potential. Finally, the argument extends 
without difficulty to higher angular momentum. 

B The complex delta functional 

In quantum mechanics, the complex delta function was originally introduced by Nakan- 
ishi [42] to describe resonances in the Lee model [43]. In mathematics, the complex 
delta function was introduced by Gelfand and Shilov [44]. The purpose of this appendix 
is to introduce the precise mathematical definition of the complex delta function and 
to show that, when the test functions are analytic, such definition coincides with the 
one given by Nakanishi. 

B.l Three definitions of the (linear) complex delta functional 

The complex delta functional has different forms depending on the properties of the test 
functions on which it acts. We shall review the three most important forms, namely 
when the complex delta functional acts on analytic functions (this form is used in this 
paper and and was introduced in [44]), when it acts on meromorphic functions (this is 
the form used by Nakanishi [42]), and when it acts on non-meromorphic functions (this 
form was introduced in [44]). When the space of test functions are analytic, as is our 
case, these three forms coincide (as they should) and can be written as in Eq. (13.151) . 
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B.l.l First definition — the test functions are analytic 



According to page 1 of Volume I of Ref . [44] , a distribution is a function that associates 
a complex number with each function belonging to a vector space: 

distribution : {Space of functions} i — > C ,„ , 

function i — > complex number . 

The functions in the "{Space of functions}" are usually called test functions. Because a 
distribution maps functions into complex numbers, they are usually called functionals. 
Such functionals can be linear or antilinear. 

A more precise definition is the following. If $ is a vector space of test functions 
endowed with a topology, a linear (antilinear) distribution F is a function from $ to C 

F ■ $ i ► C 

^ F(0) (B-2) 

such that 

(i) F is well defined, 

(ii) F is linear (antilinear), 
(m) F is continuous. 

A very important example of distribution is the (linear) Schwartz delta functional 
at a real number E. Such functional associates with each test function the value that 
takes at E: 

$E '■ $Schw 1 > C , s 

8 E {<j>) = <j>{E) , ^"^ 

where the test functions of $schw are infinitely differentiable and of polynomial falloff. 
It is straightforward to show that definition (IB. 31) satisfies the above requirements (i)- 
(iii). 

The (linear) complex delta functional is defined in a completely analogous way. 
As stated by Gelfand and Shilov [44, Vol. 2, page 85], the point E in Eq. (1B.3I) may 
be complex in the spaces of analytic functions. If $ a nai denotes a vector space of 
analytic functions at the complex point zq, then the linear complex delta functional at 
Zq is defined as a function that associates with each test function <fi the value that the 
analytic continuation of <fi takes at z : 



^zq • ^anal 



= </>{*) . (R4) 



Two important comments are in order here. First, the test functions of $ a nai must 
be analytic at z ; that is, z is not a singularity (e.g., a pole) of any 0, otherwise 
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definition (IB. 41) makes no sense. And second, the complex delta functional is completely 
specified by Eq. ( 1B.4j) because the test functions are analytic at z , and therefore one 
does not need to introduce any contour in the definition of (1B.4I) . even though one 
could use such a contour, as in Eq. ( IB. 51) below. 

Definition (1B.4|) actually fulfills the requirements («)-(m)J§ The only property that 
is conceptually challenging is (i). Because we are assuming that the test functions 
are analytic at z Q , (fi(z ) exists and is unique, which grants requirement (i). Thus, 
definition ( 1B.4j) completely, rigorously and unambiguously defines the complex delta 
functional. 



2 • - iiki u — /D C\ 

*) ( B - 5 ) 



B.1.2 Second definition — the test functions are meromorphic 

Many functions are not analytic but just meromorphic. That is, when we analytically 
continue them, they have isolated singularities ("poles") in the complex plane. At such 
poles, definition (1B.4I) makes no sense, and one has to extend it. If $ mer o is a vector 
space of meromorphic functions at Zq, the (linear) complex delta functional at z is 
defined as 

Sz '■ ^mcro 1 > C 

One can again check very easily that definition (IB. 51) satisfies requirements (i)-(iii). 
Note that because in definition (IB. 51) the test functions are meromorphic, such definition 
depends on Cauchy's theorem and on the contour used@ 

If we denote by ao the zeroth term of the Laurent expansion of <p(z) around zq, then 
definition (IB.5j) associates a with each test function 0, since 

ao = ±I dz Ji±. (B.6) 

27T1 J Z- Z 

Thus, we may write definition (IB. 51) as 

<P 1 — * 5 Z0 ((j)) = a . 



Obviously, both (IB.5j) and (IB.7j) define the same functional, because both associate 
the same complex number with the same function, even though in (IB.7I) no contour 
integral has been explicitly used. 

Now, when <f>(z) is analytic at zq, ao is simply 4>(zq). Thus, when the test functions 
are not just meromorphic but also analytic at Zq, definitions (IB. 51) and (IB. 7|) become 



3 In this paper, we omit any explicit discussion on the continuity requirement (Hi). The reason is 
that first, the continuity of the complex delta function is guaranteed by the results of [37], and second, 
continuity is not essential to our main discussion. 

4 The contour used in Eq. (|B.5j) is assumed to be a circle around zq such that the test function </> 
is analytic inside such circle except perhaps at Zq. 
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definition (1B.4j) . because in such case all these definitions associate each function with 



one and the same complex number 4>(zq). This is why, when the test functions are 
all analytic at zq, one can define the complex delta functional by way of Eq. (1B.4I) . as 
Gelfand and Shilov do in page 85, Vol. II of [44]. 

B.1.3 Third definition — the test functions are not meromorphic 



When the test functions are not meromorphic, definitions ( IB. 41) . ( IB. 51) and ( IB. 71) make 



no sense. One can still define a complex delta functional at the origin following the pre- 
scription of Gelfand and Shilov [44, Vol. I, Appendix B]. When the functions are mero- 
morphic, such definition of the complex delta functional at the origin becomes (1B.5j) 
and (1R71) . 

However, because in this paper we use test functions that are analytic at the reso- 
nance energies, we do not need to use this general definition or definition ( 1B.5I) . because 
all these definitions actually become (IB. 41) . 

B.2 Three definitions of the (antilinear) complex delta func- 
tional 

In this paper, we have used antilinear (rather than linear) functionals. We will therefore 
briefly explain how one defines such functionals for the cases considered in the previous 
section. 

The (antilinear) Schwartz delta functional at a real number E associates with each 
test function 0, the complex conjugate of the value that takes at E: 

Se ■ $schw 1 — ► C ,g x 

8 E {4>) = 0(E)* . 1 • ' 



h- 



When we write the action of 5^ as an integral operator, the kernel of such integral 
operator is Dirac's delta function: 

/■oo 

4(0) = / dE' 5{E' - E)(j>{E')* = 0(E)* . (B.9) 
Jo 

If ^anai denotes a vector space of test functions such that 0* are all analytic at 
Zo, then the antilinear complex delta functional at z Q is a function that associates with 
each test function 0, the value that the analytic continuation of 0* takes at Zq: 



®zq • ^anal 



M0)=0(^)*. (R10) 



When we write the expression for 5 Zo as an integral operator, the kernel of such integral 
operator is the complex delta function: 

POO 

lM = / dE'5(E' - z )cf>(Ey = 0( 2 *)* • (B.11) 
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When the test functions are only meromorphic and zq is one of their poles, defini- 
tion (IB.lOj) needs to be changed to 



z • =■= mero - /n i<i\ 

(B.12) 

2-20 



If we denote by ag the zeroth term of the Laurent expansion of 4>{z*)* around z , then 
definition (1B.12I) associates with each test function 0, and therefore we can write 

0~z • $mero 1 > C (B 13) 

1 — > 5 Z0 (<p) = a* . 

If the functions are not even meromorphic, we need to use the prescription of 
Gelfand and Shilov [44, Vol. I, Appendix B]. 

The same conclusions as in the previous section apply to the antilinear complex 
delta functional. When (f)(z*)* are all analytic at zq, is simply 0(zq)*. Thus, when 
the test functions are all analytic at zq, definition ( IB. 121) becomes definition (IB. 101) . 
and we are allowed to use ( IB. 101) . 



B.3 Nakanishi's definition 

Nakanishi [42] uses a slightly different version of the complex delta function. When he 
writes 6j$(<f>) as an integral operator, Nakanishi uses the following expression: 



M0) = / dE(j)(E*)*5 n (E - z K ) , (B.14) 

where 
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ME -, R) = ^(_ F i___ ; J_), (B. 15) 

and where the contour 7 is such that the integral in Eq. (1B.14I) decomposes into two 
terms. The end points of the integration paths are the same for the two terms, namely, 
and +00. The integration path for the first term, E (-]_ , passes below zr, whereas 
the integration path for the second term, E (+l_ , passes above zr. Adding the two 
terms we obtain 



/ &E<P{Ey5v{E-z K ) = ^l&E^^ = <P{z* K y. (B.16) 

J j 27T1 J E - Z R 

Thus, the distributional definition (IB. 10[) is equivalent to Nakanishi's definition flB.140 - 
(1B.16|) . because both approaches associate the same complex number, 0(z R )*, with the 



same test function, 
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